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ABSTRACT 

Variational  methods,  originally  developed  for  deriving 
stability  criteria  from  energy  principles,  are  well  suited 
for  obtaining  bounds  for  the  maximum  growth  rate  In  arbi- 
trary equilibria.   An  upper  bound  Is  derived  which  corre- 
sponds to  a  sufficient  stability  criterion,  and  a  lower 
bound  is  derived  which  corresponds  to  a  necessary  criterion 
Each  of  these  bounds  equals  the  actual  maximum  growth  rate 
in  special  cases. 
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I.   INTRODUCTION 

There  are  two  traditional  approaches  to  the  problem  of 
plasma  stability:   the  energy  principle,  and  the  normal  mode 
analysis.   The  advantage  of  the  energy  principle  is  its 
flexibility  which  allows  using  direct  variational  methods  to 
derive  stability  criteria  with  a  wide  range  of  applicability. 
On  the  other  hand,  this  flexibility  is  traded  against  infor- 
mation:  the  energy  principle  tells  us  whether  an  equilib- 
rium is  exponentially  stable  or  not,  without  giving  any  pre- 
cise information  about  the  nature  of  instabilities.   From 
the  point  of  view  of  controlled  fusion,  this  is  a  serious 
drawback  because  exponential  stability  (which  is  difficult 
to  achieve  in  practice)  is  an  unnecessarily  restrictive  re- 
quirement.  For  instance,  unstable  modes  with  small  growth 
rates  may  be  ignored  for  practical  purposes. 

A  normal  mode  analysis  yields  more  information  than  the 
energy  principle,  but  is  so  involved  that  it  has  been  car- 
ried out  only  for  special  equilibria.   As  a  consequence,  no 
results  seem  available  which  are  general  enough  to  yield  an 
Intuitive  picture  of  the  relation  between  the  geometry  of 
the  equilibrium  and  the  nature  of  instabilities. 

In  the  present  paper  an  approach  is  used  which  is  in- 
termediate between  the  energy  principle  and  the  normal  mode 
analysis  regarding  both  complexity  and  obtainable  informa- 
tion.  This  approach  is  based  on  a  modified  energy  prin- 
ciple.   Using  techniques  known  from  numerous  energy  prin— 


clple  analyses,  it  yields  bounds  for  the  maximum  growth  rate. 
More  precisely,  any  process  leading  to  a  sufficient  stabil- 
ity criterion  according  to  the  energy  principle  yields  an 
upper  bound,  and  any  process  leading  to  a  necessary  crite- 
rion yields  a  lower  bound.   The  huge  amount  of  work  already 
invested  In  evaluating  the  ordinary  energy  principle  thus 
gains  new  practical  significance. 

Even  though  this  method  applies  to  any  systems  governed 
by  a  self-adjoint  set  of  equations  (e.g.,  the  guiding  center 
plasma),  we  restrict  ourselves  to  ideal  magnetohydrodynamlcs . 
We  neither  try  to  give  an  account  of  what  one  can  achieve 
with  this  method  (which  would  Involve  reviewing  the  existing 
energy  principle  calculations),  nor  do  we  try  to  give  an- 
swers to  practical  questions  arising  from  present  fusion  ex- 
periments (which  would  Involve  restricting  consideration  to 
certain  modes  which  have  been  observed),  but  rather  give  two 
particularly  simple  (though  somewhat  academic)  illustrations, 
which,  however,  include  arbitrary  geometries. 


II.   THE  MODIFIED  ENERGY  PRINCIPLE 

The  equations  of  ideal  magnetohydrodynamlcs ,  when  lin- 
earized about  some  static  equilibrium,  are 

%   ^  =  -^Pl  Ml  ^  ^o  Mo  ^  ^1   '  (1) 

-rr  =  ^^^1  (^1  ^  ?o)    '  (2) 

3p 

W-  =  -ll   '   ^Po  -  ^Po  ^^^  ^1   '  ^3) 

where  the  quantities  with  a  subscript  one  (perturbing  quan- 
tities) are  viewed  as  unknowns,  and  those  with  a  subscript 
zero  are  known  functions  in  space,  subject  to  the  equilib- 
rium equations 

-io  ^  ^o  =  ^Po  '  io  =  ^^^1  ?o  '  div  ^^  =  0  .  (k) 
For  a  normal  mode  analysis,  one  Fourier-decomposes  the  per- 
turbing quantities  in  time,  thus  replacing  3/8t  by  iu.   Elim- 


inating ^^,    jj^,  and  p^. 


one  then  obtains 
w^Pj^  =  ^Ji   ,  (5) 

^j^  =  ^  X  curl  curl  (j^  x  ^)  -  ^  x  curl  (^   x  b) 

-  V()^-Vp  +  YP  div  ]^)  ,  (6) 
where  u  is  the  Fourier  transform  of  v-,  ,  and  subscripts  have 
been  omitted.  Since  ;^  is  a  symmetric  (i.e.,  formally  self- 
adjoint)  operator  ' -^  in  the  Hilbert  space  of  square  Integra- 
ble  vector  functions  ]^  satisfying  appropriate  boundary  con- 
ditions (e.g.,  u   =  0  at  a  toroidal  boundary),  the  spectrum 

2 
of  Eq.  (5)  is  real,  and  its  point  eigenvalues  w  are  obtain- 


able  from  the  variational  principle 


?   =  W/N  ,     6(W/N)  =  0  ,  (7) 


where 


An  equilibrium  is  called  "exponentially  unstable"  if 
the  equations  of  motion  (l)-(3)  have  a  solution  whose 
kinetic  energy,  defined  by 

K(t)  =  N[v^(^,t)] 
satisfies  the  inequality 

K(t)  >^  K  exp  (2at)  (8) 

for  some  positive  numbers  o  and  K^  and  for  large  enough  t. 

This  is  the  case  if,  and  only  if,  the  spectrum  of  Eq.  (5) 

2      2 
extends  below  the  origin.   For  instance,  if  to   =  -a   <  0  is 

a  point  eigenvalue,  equality  holds  in  (8),  and  a  is  the 

growth  rate  of  an  unstable  eigenmode.   The  "maximum  growth 

rate"  is  defined  as  the  supremum  of  all  numbers  a  for  which 

an  inequality  of  the  form  (8)  holds,  or  equlvalently ,  as 

2 
minus  the  smallest  number  w  belonging  to  the  spectrum. 

Since  the  spectrum  always  includes  the  origin,  the  maximum 

growth  rate  is  at  least  zero.   In  an  exponentially  stable 

equilibrium  it  is  zero.   It  may  also  be  infinite.   However, 

this  cannot  happen  unless  the  mass  density  is  zero  somewhere 

2 

in  the  plasma  domain,  or  the  pressure  gradient  is  unbounded. 

It  should  be  pointed  out  that  the  maximum  growth  rate  need 
not  correspond  to  an  eigenvalue.   In  axially  symmetric  equi- 
libria with  ^'J^  =  0,  for  instance,  it  corresponds  to  an  ac- 
cumulation point  of  eigenvalues  which  is  not  an  eigenvalue 
itself,  thus  belonging  to  the  continuous  spectrum. 


The  ordinary  energy  principle  asserts  that  positivity 

of  the  quadratic  functional  W  is  necessary  and  sufficient 

for  exponential  stability.   Correspondingly,  the  modified 

energy  principle  asserts  that  positivity  of  the  functional 

W  =W  +  a^N  (9) 

o 

is  necessary  and  sufficient  for  o    >   o         .        Since  this  is 
"^  —  max 

equivalent  to 

general  theorems  (viz.,  the  fact  that  the  operator  ^  has  a 
self-adjoint  extension  whenever  this  infimum  exists;  and  the 
spectral  theorem  for  self-adjoint  operators)  can  be  used  to 
prove  it.   However,  a  more  elementary  demonstration  is  more 
to  the  point  in  the  present  context. 

Necessity  of  the  modified  energy  principle  follows  from 
a  demonstration  of  necessity  of  the  ordinary  energy  prin- 
ciple :   One  introduces  a  displacement  vector  ^  by 

It  ^^^>^^  =  ^1^^'^^   ' 
so  that  Eqs.  (2)  and  (3)  can  be  integrated  to  yield 

^1  =  curl  (^  X  ^o^  +  S 

Pi  =  -^-^Po  -  ^Po  di^  ^  +  ^   ' 
where  ^(^)  and  ^i)^)    are  constants  of  integration.   Re- 
stricting consideration  to  perturbations  which  are  acces- 
sible from  the  equilibrium  through  an  actual  motion,  one 
then  puts  ^  =  0  and  P  =  0  to  obtain 

pa^^/at^  +  ^^  =  0  .  (11) 

From  this  equation,  the  energy  conservation  law 

N[8^/8t]    +   WCg    =    E   =   const  (12) 


and  the  virlal  theorem 

2-2  N[^]  -  2N[9^/8t]  +  2W[^]  =  0 
dt 

are  derived  by  taking  its  scalar  product  v;ith  8^/3t  and  with 
^j   respectively.   It  follows  from  these  two  relations  that 
a  solution  of  Eq.  (11)  exists  whose  kinetic  energy 
K(t)  =  N[3j^/3t]  satisfies  the  inequality  (8)  if  a  trial 
vector  function  j^   can  be  found  such  that  W  [jj]  <  0. 

Sufficiency  of  the  ordinary  energy  principle  is  widely 
believed  to  follow  from  the  energy  conservation  (12)  because 
W  _>  0  implies  N[v,]  <_   E,  and  hence  boundedness  of  K(t). 
However,  this  is  Incomplete  because  Eq .  (11),  and  hence  also 
the  energy  conservation  law  (12),  hold  only  for  the  subclass 
of  accessible  perturbations.   Indeed,  it  has  been  shown  that 

the  kinetic  energy  (of  inaccessible  perturbations)  can  grow 

7 
quadratically  in  time  even  if  W  _>  0.  '   Sufficiency  of  the 

ordinary  energy  principle  for  exponential  stability  to  arbi- 
trary perturbations  follows,  in  the  special  case  a  =  0,  from 
the  following  demonstration  of  sufficiency  of  the  modified 
energy  principle.   Differentiating  Eqs.  (2)  and  (3)  with  re- 
spect to  time,  and  substituting  the  resulting  expressions 
into  Eq .  (1)  one  obtains 

p9^;^/3t^  +   IX   =    0       .  (13) 

The  velocity  vector  of  arbitrary  perturbations  thus  satis- 
fies the  same  equation  as  the  displacement  vector  of  acces- 
sible perturbations.   However,  not  every  solution  of  Eq .  (13) 
corresponds  to  a  solution  of  the  original  system  (1)— (3)  be- 
cause of  the  differentiation.   Equation  (13)  thus  could  not 


be  used  for  proving  necessity  (which  involves  exhibiting  one 
solution  of  the  original  system  with  certain  properties); 
but  it  is  suitable  for  proving  sufficiency  (which  involves 
showing  that  all  solutions  have  a  certain  property).   Equa- 
tion (13)  implies  the  conservation  law 

ay 

N[x^]  +  W[v]  =  C  =  const   ,  (14) 

which  is  analogous  to  energy  conservation,  but  which  cannot 
be  identified  with  it  (for  inaccessible  perturbations,  there 
is  no  energy  conservation  law  involving  only  first  order 
quantities).   Assuming  now  that  W-  >_  0  for  some  a  >  0,  we 

have 

N[9)^/8t]  1  C  +  B^N[v]   .  (15) 

Making  use  of  the  relation  between  y;  and  ^   we  have 

11=  (;^,pa^/9t)   ,  (16) 

and  making  use  of  Schwartz'  inequality  we  have 

(;^,p8v/9t)^  <  i<N[y]N[3;//3t]   .  (17) 

Combining  the  three  relations  we  conclude  that  the  kinetic 
energy  satisfies  the  differential  inequality 

II  <  2(5^k2  +  CK)^/2   .  (18) 

-2222  - 

Since  a  K  +  CK  <  a  K  for  large  enough  K  whenever  a  >  a,  an 

inequality  of  the  form  (8)  cannot  hold  if  a  >  a.   Hence 

a    £  ^-   Fo^  0=0  the  inequality  (l8)  implies  that  K  is 

bounded  by  a  quadratic  function  of  time,  which  because  of 

the  above  mentioned  possibility  of  a  quadratic  growth  of  K, 

shows  that  the  present  estimate  is  optimal. 


Determining  a  ^   essentially  amounts  to  a  normal  mode 
max 

analysis.   The  difficulties  of  such  an  analysis  are  avoided 

if  one  merely  wants  to  calculate  bounds  for  a    .   A  lower 

max 

bound  can  be  obtained  from  a  necessary  criterion  for  W  _>  0 
by  determining  the  smallest  number  a  for  which  the  criterion 
is  satisfied.   Correspondingly,  an  upper  bound  can  be  ob- 
tained from  a  sufficient  criterion.   Criteria  for  W   >  0  are 

a  — 

obtained  in  the  same  manner  as  those  for  W  >_  0 :   Narrowing 
the  admissibility  class  of  displacements  ^  or  dropping  neg- 
ative terms  in  W  yields  a  necessary  criterion,  while  widen- 
ing the  class  or  dropping  positive  terms  yields  a  sufficient 
criterion.   In  particular,  one  can  use  various  procedures 
already  existing  in  the  literature,  thus  obtaining  bounds 
for  a    which  correspond  to  well  known  stability  criteria. 


III.   LOWER  BOUND  FOR  THE  MAXIMUM  GROWTH  RATE 

We  consider  arbitrary  toroidal  equilibria  [Eqs.  (4)], 
in  which  the  plasma  Is  either  In  contact  with  a  perfectly 
conducting  rigid  well,  or  surrounded  by  a  nonconducting  vac- 
uum which  terminates  at  a  conducting  well.   The  pressure 
surfaces  form  sets  of  nested  torolds.   Since  our  analysis 
will  turn  out  to  be  local  at  magnetic  field  lines,  we  may 
ignore  the  complications  arising  from  disconnected  pressure 
surfaces,  thus  restricting  attention  to  just  one  set  of 
nested  toroidal  surfaces. 

These  surfaces  can  be  labelled  by  their  volume  v.   The 
fluxes  \l)   and  x  the  long  and  the  short  way  of  the  magnetic 
field,  as  well  as  the  fluxes  ({>  and  n  of  the  current  density 
can  then  be  Introduced  as  functions  of  v.   These  fluxes  are 
related  to  the  pressure  through 

p  =  ri"!^  -  <^X   ,  (19) 

where  dots  denote  derivatives  with  respect  to  v.   Two  angle- 
like  variables  6  and  i.   increasing  by  unity  once  the  short 
and  the  long  way  around,  respectively,  can  be  constructed 
such  that  (v,9,c)  is  a  coordinate  system  with  Jacobian  unity, 
Vv(ve  X  Vc)  =  I,  and  such  that  the  fields  ^   and  ^  have  the 
simple  representations 

Here  ^.  =  9^^/91  d  =  v,0,^)  are  the  covarlant  basis  vectors. 
These  satisfy  the  relations 

^^  =  V0  X  V^  ,     ^g  =  VC  X  Vv  ,     ^^  =  Vv  X  V0  , 

^^  =  ^9  ^  ^C  '     ^Q  =  ^^  ^  ^v  '     ^^  =  ^v  "  ^6  • 


The  metric  tensor  g.,  =  ^-'^y.   is,  apart  from  the  fluxes,  the 

only  equilibrium  quantity  which  will  enter  our  results.   It 

can  be  determined    from  the  fields  ^   and  ^  (and  their 

fluxes)  through  the  relations  (20),  and  through 

e^  =  W/IWI^  -  k^^g  +  k^e^^  , 

where  the  coefficients  k.  are  determined  by  integrating  the 

relations 

Vv  X  Vk^  =  curl  {^^    X  Vv/lVvl^)   . 

If  all  magnetic  field  lines  are  closed,  the  rotation 

number  y  =  x/^  is  a  rational  constant.   Interchanges,  that 

is  displacement  vectors  which  leave  the  magnetic  field 

unperturbed, 

Q  =  curl  (c  X  B)  =  0   , 

may  then  be  unstable,  and  the  corresponding  necessary  sta- 
bility criterion  is  I  _>  0 ,  where 

1    =   -    n    +   ^  n^    ,  f2  =  -  ^q/q  ,     q  =  (f  d£/B  .    (21) 

Considering  only  interchanges  in  W  ,  we  now  derive  a  lower 
bound  for  a     (denoted  by  o-j.)  which  corresponds  to  the  sta- 
bility criterion  I  >^  0 . 

Since  derivatives  along  field  lines  will  play  an  impor- 
tant role,  it  is  convenient  to  replace  the  coordinate  9  by  a 
linear  combination  of  6  and  r,   which  is  constant  .on  a  field 
line.   Such  a  combination  is 

w  =  6  -  PC   . 
In  the  coordinate  system  (v,w,c),  B  =  ii^    ,    and  B'V  =  iid/'^r,. 
(Note  that  in  the  new  system  e   is  not  the  same  vector  as  in 
the  old  one,  whereas  e      remains  unchanged,  and  ^  =  ^q • ) 
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Introducing  contravariant  components  of  ^  by  writing 
one  finds  that  ^  =  0  Is  equivalent  to 

For  Interchanges,  the  first  two  components  of  ^  are  thus 

constrained  by 

X=i|i,     Y  =  -iii  (22) 

l(;  3w  '  ^  9v 

where  $  is  an  arbitrary  function  of  v  and  w,  periodic  in  w. 

Our  variational  functional  can  be  written  as 

W^  =  Wp  +  Wg  +  W^  +  a^N  ,  (23) 


where 

W 

2 


(2i|) 


I  I  d^iCl^l^  +  (^x^).g  +  (^.Vp)  div  ^ 

+  YpCdiv  ^)2]   . 

The  surface  term  W„  and  the  vacuum  term  W   (which  are  pre— 
o  V 

sent  only  if  a  vacuum  surrounds  the  plasma)  are  not  needed 
because  they  are  zero  if  one  considers,  as  we  shall  do,  dis- 
placements which  are  nonzero  only  within  the  plasma.   For 
Interchanges,  Eq .  (23)  then  becomes 

W^  =  I  I  dvdwdcClX^  +  YP(||)^  +  a^pl^l^]   .     (25) 

Minimizing  this  with  respect  to  Z  leads  to  the  Euler 

equation 

Ypa^z/ac^  =  o^pig^^X   +  g^^Y  +  g^^Z)   .        (26) 

This  is  an  ordinary  differential  equation  for  Z  with  v  and  w 

as  parameters,  subject  to  the  periodicity  condition 

Z(C+n)  =  Z(c) 

(slnglevaluedness  of  Z).   Here  n  is  the  number  of  circuits 

the  long  way  of  a  (closed)  field  line  (y  =  m/n,  m  and  n  are 


11 


\-k\ 


relative  primes).   Since  the  homogeneous  equation  (obtained 

by  putting  X  =  Y  =  0)  has  no  nontrivial  periodic  solution  if 

2 
a   >  0,  the  full  equation  has  exactly  one  such  solution  for 

every  given  X  and  Y.   Since  we  are  not  able  to  write  down 
this  solution  explicitly,  we  have  to  take  recourse  to  esti- 
mates and  limiting  cases.   However,  it  is  convenient  to  dis- 
cuss the  procedure  leading  to  the  criterion  for  W   >  0  in 
general  terms:   Since  X  and  Y  are  independent  of  c,  the  so- 
lution of  Eq.  (26)  is  a  linear  combination  of  X  and  Y  with 
coefficients  depending  only  on  equilibrium  quantities  (and 
a).   Inserting  this  solution  into  Eq.  (25),  and  performing 
the  ^-integration,  one  obtains  an  expression  of  the  form 

dvdw(W^^X^  +  2W^yXY  +  WyyY^)  ,         (27) 
where  the  components  of  the  matrix  W  „  are  functions  of  v, 
w,  and  o. 

The  posltivity  of  the  matrix  W  „,  i.e. 

W    >0,     W    >0,     WW   -W^>0 
XX  -    '      yy  -  XX  yy     xy  - 

is  necessary  and  sufficient  for  W  ^0.   Because  of  the  con- 
straint (22),  only  sufficiency  is  obvious.   To  demonstrate 
necessity,  we  observe  the  constraint  by  writing 

W   =  ^   dvdw  -^[W   l:r— 1  -2W   TT—  7—  +  W   Itt—J  J  , 
a    2  J       :2'-  xx^aw-"     xy  3w  3v    yy^9v^    ' 

and  assume  that  the  matrix  W  „  is  nonpositive  in  some  sub- 

region  of  the  (v,w)  plane.   After  performing  a  change  of 

variables  such  that  W   is  zero  at  some  point  (v^,w  )  of  this 
xy  ^       o'  o 

sub-region,  we  can  conclude  that  one  of  the  diagonal  ele- 
ments, say  W   ,  is  negative  in  a  neighborhood  of  this  point. 
A  function  $  which  is  nonzero  only  in  an  c-neighborhood  of 
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this  point,  and  whose  v-derivative  is  large  compared  to  its 

w-derlvative  (e.g.,  which  is  a  function  of  the  single  vari— 

2     2       2 
able  (v-v  )   +  e  (w-w  )  ),  will  then  render  W  negative  for 

sufficiently  small  e.   For  the  following  estimates  we  thus 

note  that  the  constraint  (22)  may  be  ignored  whenever  this 

is  convenient. 

We  now  derive  both  a  lower  bound  and  an  upper  bound  for 
Oy,    and  then  show  that  each  of  these  bounds  is  approached  in 
limiting  cases.   For  the  lower  bound,  we  further  restrict 
the  class  of  displacements  by  3Z/9c  =  0.   Equation  (25)  then 
can  be  written  as 

W^  =  I  j  dvdw  W.^X^X^ 

where  r'-=X,  X^  =  Y,  x3  =  Z, 

"l  o  o  >    2 
O  O  O  )  +  0  <pg. 
^O  O  O 

and 


=h:— =n 


61 
a^  ...  =  -  Y 

0 

is  the  usual  field  line  average.   The  positivity  of  the  3><3 


matrix  W.,  is  easily  seen  to  be  equivalent  to  the  positivity 

of  its  determinant,  and  hence  to 

I  det^^<pg^j^>  +  o^   det  <PB^^>    >  0  , 

where  det   ...  is  the  indicated  2x2  sub-determinant.   The 
wc 

smallest  value  of  a  for  which  this  holds  .  (denoted  by  a_j)    is 
a  lower  bound  for  Oj   and  is  given  by 

„       -I  det   <pg.,  > 

°?  =  ^ax    ,  .  ^^^  f    .  (28) 

-I         det  <Pgj^i^> 
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For  the  upper  bound  (denoted  by  Oj)   we  drop  the  posi- 
tive second  term  in  the  integrand  of  Eq.  (25),  thus  dealing 
with  the  same  expression  for  W  as  in  the  derivation  of  o-r.. 
But  since  there  now  is  no  constraint  on  Z,  the  condition  for 

W  is  less  restrictive,  and  leads  to 
a 

-Kpg'^Vg^  > 

0^  =  max  ^^ (29) 

det^<pg  Vg^^> 

ik 
where  g   =  Vl«Vk  is  the  contravariant  metric  tensor. 

To  show  that  the  bounds  Oj   and  o,  are  actually  ap- 
proached by  Oj,   we  solve  the  Euler  equation  (26)  in  the  two 

2  2  2 

limiting  cases  c   <<  1  and  c   >>  1,  where  c   is  a  character- 

2  2 
istic  value  of  the  quantity  n  a  pg  >/yP •   In  both  cases  the 

solution  can  be  expanded  in  powers  of  the  respective  small 
parameter.   Substituting  its  lowest  order  into  Eq .  (25),  and 
calculating  Oj   in  the  same  manner  as  before,  we  find  Oj   ^  Oj 
if  c^  <<  1,  and  a^  Z  o-^   if  c^  >>  1. 

In  terms  of  the  original  coordinates  (v,e,5),  Eqs.  (28) 
and  (29)  are 


2       -I  det   <pg.^: 
a-j-  =  ma^ 


(30) 


-2       -I<pg^^/B^>  ,^,s 

Oj   =   max  ^^p  ,  (31) 


where 


D  =  ^2(<pg^VB2><pgee/B2>  _  <pg^VB2>2] 

-2^x(<Pg^VB2><pg^W>   -    <pg^VB2><pg^W>) 
+x'(<Pg^VB2><pg^W>  -    <pg^W>2)       . 


IH 


since  a.,   <    o^   <    o        ,  a.^  is  a  useful  estimate  for  a 

-I  —  I  —  max'  -I  max 

whenever  the  criterion  I  _>  0  is  violated.   On  the  other  hand, 

o  has  no  practical  significance  in  general  because  it  can 

be  either  smaller  or  larger  than  a    .   However,  o^   ?^  a        , 

max  J.    max 

and  hence  a    <0j   if  the  magnetic  field  is  sufficiently 
close  to  the  vacuum  field  of  external  currents  (low  B).   For 
If  one  expands  the  modified  energy  principle  in  powers  of  6 
(similar  to  an  expansion  of  the  ordinary  energy  principle   ), 
one  finds,  in  the  zeroth  order,  W^°^  >  0,  and  W^°^  =  0  if, 
and  only  If,  a^°^  =  0  and  !§^°^  =  0.   Using  this  In  the  first 
order  one  exactly  repeats  the  above  calculation,  thus 
obtaining 

%ax  =  4   +  0(^')  '     4   =  °(^)   • 
The  bound  Oj   is  still  valid  if  the  magnetic  field  nas 
small  shear.   For  if  one  puts  y  =  y   +  ey-,(v),  where  y   is 
a  rational  constant,  and  then  considers  displacements  with 
Q  =  0(e),  one  finds  o-^  +  0(e)  <  a        ,  where  o.^   Is  computed 
for  the  unperturbed  closed  line  equilibrium  (e  =  0).   This 
naive  conclusion  should  not  be  Invalidated  by  the  fact  that 
the  low  shear  limit  is  singular  because  results  obtained  by 
assuming  closed  field  lines  seem  to  be  always  more  optimis- 
tic than  the  low  shear  limit  of  results  obtained  for  sheared 

12  13 
equilibria,  with  regard  to  both  the  stability  threshold   ' 

and  the  maximum  growth  rate. 

In  a  straight  circular  cylinder,  the  coordinates 

(v,e,c)  are  essentially  the  familiar  cylindrical  coordinates, 

2 
that  Isv'^VyQ'^^),    r.'^z.      In  particular,  they  are  or- 
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thogonal,  so  that  the  off-diagonal  elements  of  the  metric 
tensor  are  zero.   Since  equilibrium  quantities  are  constant 
on  field  lines  in  this  case,  the  averages  may  be  omitted  in 
Eqs.  (30)  and  (31).   Using  det  g^^^  =  det  g^^  =  1,  we  then 
find  that  Oj   =   Oj   =   o^y    and 

a^   =   max  zllllll      .  (32) 

1  P 

This  expression  agrees  with  the  maximum  growth  rate  of  res— 
onant  modes    [i.e.,  for  modes  with  k+my  =  0,  where 

^  "^   exp  l(me  +  kc)]  in  the  limit  k  ->  0°,  m  ->  «>,  if  the  latter 

2 
is  evaluated  for  p  <<  B  . 

In  general,  Oj   and  Oj   do  not  coincide,  but  are  not 
drastically  different.   In  particular,  they  both  have  the 
property  that  the  expression  to  be  maximized  vanishes  at  the 
magnetic  axis  (where  v  =  0)  because  e  'v-  1/v,  e^  '^  v,  e  %   l. 

V  Us 

From  this  we  conclude  that  violation  of  the  interchange  sta- 
bility criterion  just  near  the  axis  does  not  cause  any  seri- 
ous growth.   Since  the  mass  density  tends  to  be  small  near 
the  edge  of  the  plasma,  violation  of  the  criterion  near 
the  edge  tends  to  cause  the  largest  growth. 
These  statements  will  be  confirmed  by  the  upper  bound  for 
a    to  be  derived  in  the  next  Section. 

Let  us  finally  discuss  the  question  under  what  circum- 
stances Oj   is  close  to  Oj   or  to  a-p.   A  rough  estimate  shows 

that 

c^  ^  o^L^/^l   =  a^L^/BvJ 

where  L  is  the  length  of  a  field  line,  v^  is  the  sound  speed, 

and  V.  is  the  Alfven  speed.   Hence,  c   is  small  (and  Oj   ^  o-^) 
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if  the  Interchange  criterion  is  only  slightly  violated  or  if 
the  Alfven  speed  is  large,  and  c   is  large  (and  o^   ;^  a,)  if 
B  is  small  or  if  the  field  lines  are  very  long  (e.g.,  close 
upon  themselves  after  very  many  circuits  around  the  torus). 
For  a  family  of  slightly  unstable  equilibria  with  different  con- 
stant rotation  numbers  y,  but  otherwise  similar  parameters,  o^ 
varies  between  Oj   and  O-.   as  y  is  varied,  depending  on  whether 
y  is  a  low  order  rational  or  not.   This  is  a  rather  subtle 
geometrical  effect  which  does  not  occur  in  cylindrical  sym- 
metry.  Whether  this  effect  has  to  do  with  the  highly  com- 
plex, and  as  yet  unexplained,  behavior  of  the  Wlla  stellar— 
ator  (which  is  a  low  6,  low  shear  device  in  which  the  rota- 
tion number  can  be  varied)  remains  to  be  seen. 


17 


IV.   UPPER  BOUND  FOR  THE  MAXIMUM  GROWTH  RATE 

We  now  consider  more  general  equilibria  In  that  we 
allow  for  varying  rotation  numbers.   If  a  vacuum  surrounds 
the  plasma  we  assume,  however,  that  the  plasma  current  den- 
sity drops  to  zero  towards  the  plasma  vacuum  Interface,  so 
that  the  surface  term  In  Eq.  (23)  is  Identically  zero. 
Since  the  vacuum  term  Is  always  positive,-^  we  then  have 

+  YpCdlv  ^)2  +  a2p|^|2] 
In  all  cases. 

We  Introduce,  similar  to  those  In  Section  III,  compo- 
nents of  £  according  to 

The  Identity  -^^ 

U^4)-^  +  (^-^P)  ^i^  ^  =  -^^^  +  2(^x^^).^x  +  dlv  (xzp^) 

is  then  used  to  rewrite  Eq .  (33)  as 

+  YP(div  ^)2  +  o2pl^|2]   . 


Here 


K  =  -Q  -   '|,[x^^l^  ,     ^   =    f)i>  -    ^x 


The  quantity  Q.   reduces  to  that  In  Eq .  (21)  if  all  field 
lines  are  closed.   Dropping  the  first  and  the  third  term  in 
Eq.  (3^^),  and  using  Schwartz'  inequality 

UI^|Vv|2  >  (^.Vv)2 
for  the  last  term,  we  see  that 

K  +  o^p/|Vv|^  21  0 
is  sufficient  for  W  _^  0 .   The  expression 
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2  -KIVVI^  f^r-^ 

o     =  max  — '   I  (35) 

K  p 

2 
thus  Is  an  upper  bound  for  a    which  corresponds  to  the 

sufficient  stability  condition   K  >_  0 . 

The  stability  condition  K  >_  0  is  violated  near  a  mag- 
netic axis  (Vv=  0)  whenever  J  7^  0  there,  because  then 
Ij^x^  I   =  0(l/v).   The  bound  (35)  shows  that  the  growth  rate 


annihilates  tte  singularity  in  the  stability  criterion. 
Moreover,  the  expression  to  be  maximized  in  Eq .  (35)  van- 
ishes at  the  magnetic  axis  whenever  f^  =  0  there,  so  that  vi- 
olation of  the  criterion  K  _>  0  just  near  the  axis  It^.ds  only 
to  small  growth  rates.   This  is,  for  instance,  the  case  in 
stellarator  equilibria  (which  are  usually  defined  as  to  have 
no  toroidal  net  current,  that  is,  by  cf)  =  0).   This  indicates 
that  estimates  of  a  critical  B  in  stellarators  which  were 
obtained  by  an  evaluation  of  the  criterion  K  _>  0   at  the 
magnetic  axis,  can  be  significantly  improved  if  one  allows 
for  instabilities  with  a  small  growth  rate,  that  is,  if  one 
replaces  this  criterion  by  a^.  _<  a  ,  where  a      Is  some  given 
number  corresponding  to  practical  requirements.   More  gener- 
ally, results  arising  from  evaluating  any  stability  criteri- 
on at  the  magnetic  axis  are  subject  to  some  doubt  as  to 
their  relevance  because  one  has  to  consider  the  whole  plasma 
domain  in  order  to  obtain  an  upper  bound  for  a        . 
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To  give  a  specific  (though  trivial)  example,  we  consid- 
er the  straight  Z-pinch.   In  cylindrical  coordinates  (r,(i),z), 
the  magnetic  field  has  only  a  <}>  component  B(r),  and  the 
equilibrium  equation  is 

dr  2  dr 

Our  bounds  for  a    turn  out  to  be  given  by 

P  ,  ,,,  -B^(ft2)(f.2s)  (3g) 

2 
aj  =  max  =^^   ,  (37) 

where 

2    dr  »     ^     2 
It  has  been  shown  that  maximum  growth  appears  in  the  limit 
k  -*-  0°,  where  ^  ^   exp  i  (m(})+kz)  ,and  that  a    can  be  computed 
from 

(m^  +  X^+2f)(m^g+X^)  +  i^gX^  =  0  ,     X^  =  o^pr^/B^ 

2 
by  solving  for  a  ,  and  then  maximizing  with  respect  to  both 

m  and  r.   Maximum  growth  can  be  either  resonant  (m=0)  or 
off-resonant  (m/0),  depending  (in  an  extremely  complex  man- 
ner) on  the  functions f(r)  and  g(r).   A  sufficient  condition 

14 
for  resonant  maximum  growth  is 

f  +  2g(l+g)  <  0   ,  (38) 

and  o  is  then  given  by 


L.  -   »a.  --^^ 


Obviously,  a^^^   ;^  o-j.  if  |  f  |  <<  1  and  g  <<  1,  while  a^^^  ^  a^ 
if  g  <<  |f|  (at  the  maximizing  value  of  r) .   The  former  is 
the  case  at  low  B  (as  we  have  already  seen  in  more  general- 
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Ity),  and  the  latter  Is  the  case  if  the  pressure  gradient  i- 
large  at  the  plasma  edge. 

For  Instance,  the  pressure  profile 

P  =  Pq^^  -   —^    >  0  <  r  <  a 

a 
leads  to 

„2    ^"Po'-" 


and  f  =  -(n+2)/2  =  const,  whereas  g(r)  decreases  from  g  =  1 

at  the  axis  (r=0)  to  g  =  0  at  the  wall  (r=a)..   For  n  >_  6  the 

inequality  (38)  holds  so  that  o^^^   is  given  by  Eq.  (39). 

max 

2   2 
Since  B  /r  is  increasing,  while  f  is  a  negative  constant  and 

g  is  decreasing,  maximum  growth  appears  at  the  wall  (at 

least  as  long  as  p  is  not  Increasing),  and  a  =  a„  in  this 

max    IV 
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